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Abstract 

A remarkable thermodynamic fermion-boson symmetry is found for the canonical ensem- 
ble of ideal quantum gases in harmonic oscillator potentials of odd dimensions. The bosonic 
partition function is related to the fermionic one extended to negative temperatures, and vice 
versa. 
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1 Introduction and summary 



Finite quantum gases of identical fermions or bosons contained in harmonic os- 
cillator potentials are currently of great interest. The recent motivation stems from 
the investigation of alkali atoms (bosons) in magnetic traps and the observation of 
Bose-Einstein condensation, see e.g. [1-3]. On the fermionic side nuclei and their 
low temperature behaviour are investigated, for instance in [4,5]. If these systems 
are either only weekly interacting or close to the ground state, they may be consid- 
ered as an ideal gas of (quasi-) particles confined by a harmonic oscillator potential. 

For these finite ideal quantum gases simple expressions, e.g. for the partition func- 
tion of the canonical ensemble, could not be found up to now, but recursion relations 
could be derived [6] and used successfully [7-1 1]. 

In this article we report that a closer inspection of the canonical partition func- 
tion uncovers a surprising symmetry property which connects fermions and bosons. 
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Fermionic and bosonic iV-body partition functions are connected via analytic con- 
tinuation to negative temperatures. Thermodynamic mean values for the bosonic 
systems can be easily evaluated knowing the fermionic results and vice versa. Es- 
pecially for fermions this procedure removes the well-known sign problem because 
the bosonic partition function does not suffer from alternating signs. 



2 Symmetries of the partition function 



We consider N identical particles subject to a single-particle Hamiltonian in the 
canonical ensemble. The partition function Z N can be recursively built starting 
with the single-particle partition function [6] 



1 N 1 

Z n(P) = jyE ( ±:L ) n+1 Z ^P) Z N-n(P) , Z Q {0) = 1,P= — , (1) 

where the upper sign in the sum stands for bosons and the lower sign for fermions. 
For the case of a three-dimensional harmonic oscillator potential the single-particle 
partition function is 



[ exp(-/3f) 
m) = [l-exp(-/?M 

Introducing the abbreviation 



(2) 



y = exp(—(3fiw) 
we redefine the partition functions to depend on y, e.g. 
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The recursion formula (1) then reads 
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The partition functions Z^ are rational functions of y and can be split into factors 
[7,8] 
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2%(y) = N f / P±(y) (6) 

with Pjy(y) being polynomials in y. Obviously this leads to a recursion formula for 
the polynomials 



PUv) = E (±i)" +1 ^^rr^ ^-nd/) > ^ 

where Pq(v) = 1 and P^(y) = 1. This is already a useful result since for numeri- 
cal evaluation recursion formula (7) behaves better. 

A closer look at these polynomials opens new insight into fermion-boson symme- 
try. Consider for example N = 3: 



P+{y) = 1 + 3y 2 + 7y 3 + 6y 4 + 6y 5 + 10y 6 + 3y 7 (8) 
P 3 ~ (y) =y 2 (3 + 10y + Qy 2 + Qy 3 + 7y A + 3y 5 + ly 7 ) . 

It seems surprising that the coefficients appear in reverse order comparing one poly- 
nomial with the other. This is not mere coincidence, but the following theorem 
holds 



PZ(v) = V* N < N - 1) P£(jJ , (9) 

which can be proven by induction or with the help of an explicit representation 
as done below. The property of the bosonic and fermionic polynomials leads to a 
relation between fermionic and bosonic partition functions which is 



z + M = {-i) N z* [fj • do) 

Using the usual dependence on the inverse temperature (3 this property expresses 
itself as 



Z%{fi) = (-l) N Zj f {-($), (11) 

where the partition function with the negative argument has to be understood as the 
analytic continuation into the region of y — exp(—(3huj) > 1. In thermodynamic 
mean values like mean energy or specific heat this symmetry also shows up 



E+((3) = -E N (-(3) , C+(P) = C N (-(3). (12) 
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A straight forward application of the above result is to calculate fermionic partition 
functions and mean values by evaluating the respective bosonic ones at negative 
temperatures and thereby to avoid the sign problem. 

More generally this property is related to the fact that the single-particle partition 
function has an analytic continuation to the whole /3-axis where it is an odd func- 
tion, Z x {-$) = Then Z% can be analytically extended in the same way 
and it satisfies eq. (11). We prove this relation rewriting the explicit representation 
of Zx(P) as it is given in eq. (8) of ref. [11] in the form 



Z± = {±1) N N f} E (±1)fc+i ntoM(^ + i-n^) ; (13) 

k=0 0<m<-<n k <N llj=l n j 

where we understand n = and n k+1 = N in the upper product. If f3 is replaced 
by —f3 we obtain (k + 1) minus signs in the upper product, which either cancel the 
factor (— l) h+1 in the fermionic case or introduce it in the bosonic one and thereby 
transform the fermionic partition function into the bosonic one or vice versa. 

This theorem shows that the fermion-boson symmetry depends only on the oddness 
of Z\ and not on the form of the single-particle Hamiltonian. But the harmonic os- 
cillator potentials for odd space dimensions are the only examples we know, where 
this condition holds. 



3 Properties of the polynomials 



In order to study the thermodynamics of ideal quantum gases in harmonic oscillator 
potentials one would like to know more about the polynomials defined by (6). 
Using the explicit representation of the partition function (13) one derives 



N 



= (±1)" IK 1 "I/ 1 ) 3 (14) 
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N-l 



x E E (±i) fe+1 



o< ni < <n k < N n -=i nj nio (i - y n ^- n r 

Note, that P^, although looking like a rational function, is actually a polynomial 
since n-=o (1 - y n ^- n f divides Uf=i(l - y l f- 

In the following we have to distinguish between properties of P^, which can be 
mathematically proven and those which are only established by numerical evi- 
dence, i.e. usually checked with MATHEMATICA for several N. For simplicity 
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we are considering only bosonic polynomials now, the respective fermionic ones 
can be obtained by transformation. 

The following belongs to the second class (numerical evidence). The coefficients 
of the polynomials 



L(N) 

PUy)=T,p { n N) y n (is) 

n=0 

are non-negative and, for n > and iV > 4, monotonically increasing until a 
maximum is reached and then monotonically decreasing. The degree L(N) satisfies 



L(N) < i (3iV 2 - 7N + 10) . (16) 

The first N+ 1 coefficients p^ are independent of N. They form an approximately 
exponentially increasing sequence 



1, 0, 3, 7, 18, 39, 99, 213, 492, 1056, .... (17) 

For the independence of the first N/ 2 coefficients we know a mathematically exact 
proof, for the others we have numerical evidence. 

Several identities hold for the polynomials, e.g. 



P N (l) = (N\f (18) 
P' N (l) = ^N(N-l)(N\f 

P£(l) = —(27N 2 - 23iV - 8)N{N - l)(iV!) 2 
48 

P^'(l) = ^(27A^ 4 - 42A^ 3 - 9A^ 2 + 16A^ - 20)iV(A^ - 1)(A^!) 2 , 
where the primes denote derivatives. We sketch a proof for the first relation, 



L(N) 

Pn(1) = E P?° = (^O 2 , (19) 

n=0 

the others may be proven in a similar way. All terms with n > 1 in recursion 
formula (7) contain at least one factor (1 — y j ) 3 and thus vanish for y = 1. Hence 
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1 (1 -v N ) 3 

P »0) = - T \. \{ Pn-i(v) , (20) 



N {l-yf 

= ^(i + y + y 2 + --- + y N - 1 fPN-i(y) 



y=l 



= N 2 P N . 1 {1) 

The relation (19) then follows by -Pi(l) = 1 and induction. 
Introducing a quasi-statistical notation 



i 

« /W )) : = 77ml E P$Pf(n) (21) 



we can summarize these results as 



fM:=((n)) = ^N(N-l) (22) 

( x 2 :=((n 2 ))- /U 2 = l(2iV + 5)iV(iV-l) 

^:=(((n-^)) = ~N(N-l). 

The last result is remarkable since it implies that the "skewness" of the distribution 
of the pW, defined by [12] 



7 = ^ (23) 

is asymptotically proportional to N~ 5 / 2 and thus vanishes for N — > oo. This is 
compatible with the observation that for large AT the distribution of coefficients can 
be approximated by a symmetric bell-shaped function. 
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